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Exhaust  gas  travels  through  Cylinder  head  and  exhaust  manifold,  a  turbocharger  to  increase  engine  power, 
a  catalytic  converter  to  reduce  air  pollution  and  noise  reducer  such  as  muffler.  Inside  the  manifolds  the  exhaust  gases 
coming  out  of  multiple  cylinders  are  collected  and  then  directed  through  towards  catalytic  converters.  Catalytic  converters 
are  generally  fitted  in  series  with  the  exhaust  pipe  of  almost  all  vehicles.  Inside  the  Catalytic  converter  Harmful  gases  like 
Co,  No  etc.  are  converted  into  C02,  N2  etc.  respectively.  They  oxidize  carbon  monoxide  (CO)  to  carbon  dioxide  (C02), 
and  reduce  nitrogen  oxides  (NOx)  into  nitrogen  and  water  vapors.  This  work  is  based  on  combination  of  experimental  and 
computational  investigations  to  optimize  the  design  of  exhaust  system  from  emission  point  of  view.  The  aim  of  this  work  is 
to  optimize  design  of  exhaust  manifold  in  order  to  improve  filtration  efficiency  of  catalytic  converter  with  MC-2 
configuration  using  Cu  as  a  catalyst.  Amongst  various  constituents  of  exhaust  gas,  CO  and  NO  are  most  dangerous 
pollutants  from  Global  warming  point  of  view.  Thus  this  particular  work  is  concentrated  on  analyzing  filtration  efficiency 
for  these  two  pollutants. 

KEYWORDS:  CFX,  CATIA  V5,  Gambit  2.2,  DETCHEM,  Fluent  6.2,  SBCE,  SBCER,  LBCE,  LBCER,  SBSE,  SBSER, 
LBSE,  LBSER,  Relative  A:  F  Ratio  (X) 

INTRODUCTION 

Conversion  efficiency  of  catalytic  converter  is  the  measure  of  amount  of  CO  oxidized  to  C02  or  Amount  of  NOx 
reduced  to  N2. 

Kulal  et.al  (2012)  analyzed  three  different  configurations  of  catalytic  converters  viz.  MC-1,  MC-2  and  MC-3  for 
flow  characters  like  velocity  distribution  &  pressure  drop  to  optimize  the  design  for  minimum  particulate  matter  in  the 
exhaust.  Kulal  et.al  (2013)  analyzed  by  CFD  analysis  effect  of  manifold  geometry  on  emissions.  Kulal  et.al  (2014) 
analyzed  performance  of  catalytic  converter  with  MC-1,  MC-2  &  MC-3  configuration  with  metallic  copper  as  catalyst  and 
proved  that  catalytic  converters  of  MC-2  configuration  gives  best  filtration  efficiency  compared  to  others  for  harmful 
pollutants  like  NOxand  CO.  Kulal  et.al  (2014)  analyzed  performance  of  catalytic  converter  with  MC-2  configuration  with 
various  metal  catalysts  and  proved  that  copper  can  prove  to  be  an  effective  long  term  replacement  to  precious  &  rare  earth 
metals  like  platinum,  palladium  and  rhodium. 

This  particular  work  is  extended  discussion  of  these  4  previously  presented  research  papers  presented  by  us 
(Refer:  Reference:  13,  15,  16  &  19). 
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METHODOLOGY 

This  paper  deals  with  the  study  of  fluid  flow  inside  the  entire  exhaust  system  starting  through  exhaust  manifold; 
then  through  catalytic  converter  and  the  study  of  temperature  distribution  and  chemical  reaction  in  catalytic  converter. 
CATIA  V5  R-15  was  used  for  geometric  modeling  of  catalytic  converter.  Discretization  was  done  with  Gambit  2.2  and 
flow  was  analyzed  with  fluent.  Detchem  was  used  for  chemical  analysis. 

In  comparison  to  the  Navier-Stokes  equations,  the  diffusion  terms  along  the  direction  of  convection  are  neglected. 
This  model  has  been  proven  useful  for  higher  flow  velocities,  as  they  are  expected  in  this  application. 

We  have  considered  two  chemical  reactions.  The  oxidation  reaction  of  CO  and  hydrocarbons  and  reduction 
reaction  of  NO  were  considered.  The  hydrocarbons  were  represented  by  propylene,  which  is  easily  oxidized  hydrocarbon, 
constitute  about  80%  of  the  total  hydrocarbons  found  in  the  typical  exhaust  gas.  The  chemical  reactions  are: 

•  Oxidation  of  Carbon  Monoxide 

2CO+02  ->  2C02  (1) 

•  Oxidation  of  Hydro  Carbon 

2C3H6+902  ->  3C02+6H20  (2) 

•  Reduction  of  Nitrous  Oxide 

CO+NO  ->  C02  +  Vi  N2  (3) 
MODEL  DESCRIPTION 

We  have  considered  8  different  models  of  exhaust  manifold  for  this  particular  research. 
The  eight  models  considered  for  this  work  are 

•  Short  Bend  Center  Exit  (SBCE) 

•  Short  Bend  Side  Exit  (SBSE) 

•  Long  Bend  Center  Exit  (LBCE) 

•  Long  Bend  Side  Exit  (LBSE) 

•  Short  Bend  Center  Exit  with  Reducer  (SBCER) 

•  Short  Bend  Side  Exit  with  Reducer  (SB SER) 

•  Long  Bend  Center  Exit  with  Reducer  (LBCER) 

•  Long  Bend  Side  Exit  with  Reducer  (LBSER) 

The  Header  length  for  all  manifold  was  kept  same  i.e.  335  mm.  The  longer  bend  models  had  bend  radius  of 
100  mm  whereas  shorter  bend  models  had  bend  radius  of  48  mm. 
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Figure  1:  Catalytic  Converter 

For  the  simulations  a  commercially  available  three-way  catalytic  converter  with  a  rectangular  cross-section  was 
assumed.  The  specifications  of  the  substrate,  the  coating  and  the  canning  are  listed  in  following  tables: 

Table  1:  Model  Description 


Property 

1st 

Compartment 

2nd 

Compartment 

Cell  Density 
(Cells/  Inch2) 

144 

324 

Hydraulic 
Diameter  (mm) 

0.51 

0.41 

Uncoated  Wall 
Thickness  (mm) 

8 

6 

Wash  coat 
Thickness  (um) 

20 

20 

Porosity  (%) 

57.7 

50.8 

Substrate:  Knitted  Steel  Wire 

Coating:  Copper  (50gm/ft3) 

Inlet  Dimensions:  38  X  56  mm 

Cross  Section:  176  X  1 10  mm 

Length  of  One  Section:  153  mm 


Figure  2:  3-D  View  of  Catalytic  Converter  Figure  3:  Side  View  of  Grid 
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Figure  4:  3-D  View  of  Grid       Figure  5:  Top  View  of  Grid 


Figure  6:  3-D  View  of  Two  Compartments     Figure  7:  Assembled  Compartment 


Figure  8:  Path  of  Exhaust  Gas  through  Catalytic  Converter 
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Figure  9:  Meshed  Model 

The  eight  models  are  as  follows: 


Figure  10:  SBCE  &  SBCER 


Figure  11:  LBCE  &  LBCER 


Figure  12:  SBCE  &  SBCER 
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Figure  13:  LBCE  &  LBCER 
MATHEMATICAL  AND  NUMERICAL  MODEL 

Numerical  Model  for  the  simulation  was  developed  using  below  equations.  Flow  was  transient  and  unsteady. 
The  equations  to  be  solved  are: 
•     Conservation  of  Total  Mass 

d(rpu)  d{rpv) 


+ 


dr 


=  0 


(i) 


Conservation  of  Species  Mass 

d(rpu7l)i  d(rpvYi)_  d 
dz  dr  5r 


•     Conservation  of  Axial  Momentum 

dirQua)    dirpxni)       dp     d  (     du s 
+  — — i  =  — jr-^-  +  —  j  \.u~— 


(2) 


dz 


dz  m 


dr 


(3) 


Conservation  of  Enthalpy 


ru 


dz  dr 


Xr 


3T\  d 
6r  )  dr 


(4) 


Despite  the  actual  shape  of  the  channels,  a  cylindrical  model  is  used  in  order  to  reduce  the  numerical  complexity. 

Catalytic  reactions  at  the  surface  are  taken  into  account  in  terms  of  the  diffusive  mean  flux  jt  at  the  gas  phase  wash 
coat  interface  in  Equation  2. 


(5) 
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*i = jy*kfk  n  vft 
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The  boundary-layer  equations  can  be  solved  in  a  single  sweep  of  integration  along  the  axial  direction  by  a 
method-of-lines  procedure  using  an  adaptive  integration  step  size.  The  radial  derivatives  are  discretized  using  a 
finite-volume  method.  The  resulting  differential-algebraic  equation  system  is  integrated  using  the  semi-implicit 
extrapolation  solver  LIMEX. 

Using  this  steady-state  channel  model,  a  transient  simulation  of  the  thermal  behavior  of  the  entire  catalytic 
converter  is  performed  by  DETCHEMMOnolith-  The  solid's  temperature  field  is  described  by  a  three-dimensional 
conservation  equation: 


dt 


(9) 


q  =  -a  ■  2nr  A 


or 


(10) 


Boundary  Conditions 


The  simulation  and  analysis  was  carried  out  at  different  exhaust  temperatures  ranging  between  250  degree  Celsius 
&  1000  degree  Celsius.  Also  we  considered  3 -different  A:  F  ratios: 

Table  2:  A:  F  Ratio 


Rich  Mixture 

Relative  A:F  Ratio  (lQx)  0.5 

Nearly  Stochiometric  Mixture 

Relative  A:F  ratio  (Xox)  0.9 

Lean  Mixture 

Relative  A:F  ratio  (kox  1.5 

Table  3:  Applied  Boundary  Conditions 


Entity 

Zone 

Zone  Type 

mass  flow  inlet 

boundary 

Inlet 

pressure  outlet 

boundary 

Outlet 

boundary  1 

boundary 

wall-reduction 

boundary  2 

boundary 

wall-oxidation 
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The  exhaust  gas  mass  composition  was  determined  experimentally  for  3  A:  F  ratios.  The  engine  specifications 
were  as  follows: 

Table  4:  Engine  Specifications 


Engine 

4  Stroke  4  Cylinder  SI 
Engine 

Make 

Maruti-Suzuki  Wagon-R 

Calorific  Value  of  Fuel  (Gasoline) 

45208  KJ/Kg-K 

Specific  Gravity  of  Fuel 

0.7 

Bore  and  Stroke 

69.05mm  x  73.40mm 

Swept  Volume 

HOOcc 

Compression  Ratio 

7.2:1 

Dynamometer  Constant 

2000 

Diameter  of  Orifice 

29mm 

Coefficient  of  Discharge 

0.65 

The  average  Composition  of  exhaust  gas  was  found  out  experimentally  at  various  air  fuel  ratios  by  flue  gas 
analysis.  The  air  fuel  ratio  was  varied  and  the  proportion  of  various  constituents  of  exhaust  gas  was  calculated  from  that. 

Table  5:  Exhaust  Gas  Composition 


%  CO 

%  NOx 

%  co2 

%  H20 

%  HC 

%N2 

%  o2 

Others 

Total 

A,  =1.5 

0.01484 

0.00170 

8.406 

18.126 

0.754 

71.814 

0.844 

0.03945 

100 

X  =0.9 

0.01534 

0.00167 

7.762 

18.158 

0.776 

71.986 

0.87 

0.43098 

100 

X  =0.5 

0.01588 

0.00155 

7.808 

18.078 

0.926 

71.924 

0.79 

0.45656 

100 

RESULTS 

The  exhaust  gas  composition  provided  composition  of  exhaust  gas  at  the  entry  point  of  the  model  Table  5.  Each  of 
the  eight  models  as  discussed  earlier  consisted  of  4  inlets  and  one  outlet.  The  chemical  reaction  and  flow  taking  place 
inside  the  entire  exhaust  system  was  simulated  as  discussed  earlier.  The  outlet  composition  was  evaluated  a  result  of  this 
simulation. 

Once  this  result  was  obtained,  the  conversion  efficiency  of  CO  &  NO  was  evaluated  for  each  model. 


Conversion  efficiency  of  "X"  = 


?     Inlet  mass  fract  ion  of  "X'  -  On  tlst  mass  fi-act  ion  of  "X" 


fnlst  Mass  Frac tianofX 


xlOO 


Thus  similarly  we  carried  out  24  simulations  on  each  model  by  considering  8  simulations  at  8  different  exhaust 
temperature  between  300  to  1000  °C  for  these  simulations  {At  each  of  these  eight  temperatures;  3  simulations  for 
3  different  A:F  ratios  Table  2} 

Thus,  after  evaluating  such  conversion  efficiencies  for  all  simulations  on  all  models  which  resulted  in  series  of 
values;  these  obtained  values  were  plotted  on  graphs  shown  below  Taking  1  A:F  ratio  at  time  and  plotting  conversion 
efficiency  VS  exhaust  temperature  for  all  models  resulting  in  below  6  plots  (3  for  3  A:F  ratios  and  such  two  sets  for  CO 
and  NOx  respectively). 
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Figure  14:  CO  Conversion  Efficiency  vs  Temperature  for  Different  Models  (k  =0.9) 
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Figure  15:  NOx  Conversion  Efficiency  vs  Temperature  for  Different  Models  (k  =0.9) 


-convers 

con  vers 
con  vers 
con  vers 
con  vers 
con  vers 
convers 
-convcrs 


ion  effici 
ion  effici 
ion  effici 
ion  effici 
ion  effici 
ion  effici 
ion  effici 
ion  effici 


cncy  (SBCE) 
ency  (SBCER) 
ency  (LBCE) 
ency  (LBCER) 
ency  (SBSE) 
ency  (SBSER) 
ency  (LBSE) 
ency  (LBSER) 


o  ?on        400         r.oo         roo        iooo  i?oo 


Figure  16:  CO  Conversion  Efficiency  vs  Temperature  for  Different  Models  (k  =1.5) 
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Figure  17:  NOx  Conversion  Efficiency  vs  Temperature  for  Different  Models  (k  =1.5) 
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Figure  18:  CO  Conversion  Efficiency  vs  Temperature  for  Different  Models  (k  =0.5) 
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Figure  19:  NOx  Conversion  Efficiency  vs  Temperature  for  Different  Models  (k  =0.5) 
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CONCLUSIONS 

From  Figure:  14-19  it  can  be  easily  concluded  that  conversion  Efficiency  for  NOx  was  found  to  be  highest  with 
LBCE  exhaust  manifold.  Whereas  conversion  Efficiency  for  CO  was  found  to  be  highest  with  LBCE  &  SBCE  exhaust 
manifold. 

K.S.  Umesh,  V.K.  Pravin  &  K.  Rajagopal  in  their  "CFD  Analysis  of  Exhaust  Manifold  of  Multi-cylinder  SI 
Engine  to  Determine  Optimal  Geometry  for  Reducing  Emissions"  (Refer:  References:  16)  predicted  that  LBCE  is  indeed  a 
best  model  from  emission  point  of  view  by  analysis  of  exhaust  velocity  and  back  pressure  contours  obtained  as  a  result  of 
flow  simulations  through  the  same  eight  models  under  discussion. 

Thus  our  work  justifies  those  claims  by  confirming  the  results.  Thus  we  conclude  that  LBCE  is  indeed  best  design 
from  emission  point  of  view  amongst  models  under  discussion 

Also  another  fact  that  was  observed  from  Figure  14-19: 

4  4 
X  Area  under  co  nversi  oi  efficiency  \  Ars  a  under  conversion  efficiency 

/  J  curve  for  a,  I  basEc  models  (i.e.  with  reduced      /  1  cutve  for  al.  modEfEed  rrodels  (E,  e,  wEtho  it  reducer) 

Please  note  that  we  are  referring  SBCE.  SBCE,  LBCE  and  LBSE  as  basic  models  and  same  4  models  attached 
with  reducer  as 

Modified  models  (i.e.  SBCER,  SBSER,  LBCER,  LBSER). 

Thus  we  can  conclude  that  addition  of  reducer  reduces  conversion  efficiency  for  models  under  discussion. 

Above  stated  conclusions  and  Figure  14-19  indicate  that  conversion  efficiency  of  the  various  pollutants  is 
dependent  upon  symmetry  of  the  design  in  some  or  the  other  way. 

Thus  to  find  out  the  exact  relationship  between  the  symmetry  of  the  design  and  conversion  efficiency  of  various 
we  have  evaluated  design  symmetry  of  each  design  in  terms  of  mean  deviation  of  the  lengths  of  path  travelled  by 
exhaust  gas  stream  coming  out  of  each  cylinder  up  to  outlet  of  exhaust  manifold 

For  Example 

Consider  a  SBCE  model  shown  below: 
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Y.'J 


X 

G. 


Met-1 

Met-2 

Mefc-3 

IdLet-4- 

Figure  20:  SBCE  Exhaust  Manifold  Block  Diagram 


Let, 


Li  =  Path  travelled  by  exhaust  gas  from  inlet  "i"  to  Outlet 
Thus  for  above  SBCE  model  under  consideration 
Ll=  X+  Y  +  (Y/2)  +  Z 

(Refer  figure  27  for  further  information  regarding  these  Dimensions) 
LI  =  L3=  {(7)  x  48}  +  85  +  42.5  +  220  (Dimensions  are  in  mm) 

L2  =  L4=  {0  x  48}  +  42.5  +  220  (Dimensions  are  in  mm) 
Lavg  =  " 


ALi  =  \l(Li  -  Lavg)  (Where  "i"  =  1  to  4) 


AL: 


This  "AL"  can  be  treated  as  measurement  of  symmetry  of  exhaust  manifold.  We  have  evaluated  design  symmetry 
of  all  the  models  under  consideration  as  follows: 

Table  6:  Evaluation  of  Design  Symmetry  (AL) 


Model 

LI 

L2 

L3 

L4 

Lavg 

AL 

SBCE 

422.86 

337.86 

337.86 

422.86 

380.36 

0.223472 

SBCER 

422.86 

337.86 

337.86 

422.86 

380.36 

0.223472 

LBCE 

504.5 

419.5 

341 

504.5 

442.375 

0.307624 

LBCER 

504.5 

419.5 

341 

504.5 

442.375 

0.307624 

SBSE 

337.86 

337.86 

385.18 

507.86 

392.19 

0.354523 

SBSER 

337.86 

337.86 

385.18 

507.86 

392.19 

0.354523 

LBSE 

419.5 

419.5 

504.5 

589.5 

483.25 

0.291685 

LBSER 

419.5 

419.5 

504.5 

589.5 

483.25 

0.291685 
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One  can  easily  observe  that  all  the  basic  and  corresponding  modified  models  have  same  design  symmetry  factor. 
Thus  we  have  evaluated  average  conversion  efficiency  of  those  two  models  in  below  plots  in  which  we  have  plotted 
variation  of  conversion  efficiency  with  AL  at  various  temperatures  and  A:F  ratios.  (One  can  note  that  these  plots  indicate 
same  variation  indicated  by  above  plots;  but  instead  on  8  curves  we  would  be  having  only  4  curves  per  graph) 

Sum  of  C  onv  erst  on  eff  ici  eitcy  r'o  r  mad  el  s  with  ^11  =X" 


Conversion  efficiency  for  models  with  "AL  =X"  = 
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Figure  21:  CO  Conversion  Efficiency  vs  Temperature  for  AL  (k  =0.9) 


Figure  22:  NOx  Conversion  Efficiency  vs  Temperature  for  AL  (k  =0.9) 
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Figure  23:  CO  Conversion  Efficiency  vs  Temperature  for  AL  (k  =1.5) 


Figure  24:  NOx  Conversion  Efficiency  vs  Temperature  for  AL  (k  =1.5) 


Figure  25:  CO  Conversion  Efficiency  vs  Temperature  for  AL  (k  =0.5) 
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Figure  26:  NOx  Conversion  Efficiency  vs  Temperature  for  AL  (k  =0.5) 

From  above  plots  (Figure  21-26)  we  can  conclude  that  both  CO  Conversion  Efficiency  and  NOx  Conversion 
efficiency  varies  inversely  as  AL. 

Hence  the  inference  is;  for  higher  conversion  efficiency  for  both  Co  and  NOx  the  value  of  AL  must  be  close 
to  zero  as  far  as  possible. 
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NOTATIONS 

•  a  Concentration 

•  cp  Specific  Heat 

•  fk  Dependency  Function 

•  hi  Enthalpy  of  Species  i 

•  h  Enthalpy  of  the  Mixture 

•  p  Pressure 

•  q  Heat  Source  term 

•  ji  Diffusive  Flux  Including  Surface  Flux 

•  kfk  Reaction  Rate 

•  r  Radial  Spatial  Coordinate 

•  A  Face  Area 

•  Ak  Pre-Exponential  Factor 

•  Eak  Activation  Energy 

•  Fcat/geo  Ratio  of  Catalytic  to  Geometric  Surface  Area 

•  Ks  Number  of  Surface  Reactions 

•  Mi  Molar  Mass 

•  Ng  Number  of  Gas-Phase  Species 

•  Ns  Number  of  Surface  Species 

•  R  Gas  Constant 

•  T  Temperature 

•  Aox  Relative  A:F  Ratio 

•  T  Time 

•  u  Axial  Velocity 

•  v  Radial  Velocity 

•  w  Normal  Velocity 

•  w  Mean  Normal  Velocity 

•  Yi  Mass  Fraction  of  Species  i 
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z 

Axial  Spatial  Coordinate 

O  7 

Temperature  Exponent 

zik 

Coverage  Dependent  Change  of  Activation  Energy 

Wash  Coat  Effectiveness  Factor 

X 

Thermal  Conductivity 

y 

Uniformity  Index 

fi 

Viscosity 

[iik 

Coverage  Dependent  Change  of  Reaction  Order 

vik 

Stoichiometric  Coefficients 

P 

Density 

a 

Channel  Density 

0/ 

Surface  Coverage 
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Figure  28:  SBCER 
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Figure  31:  SBSE 
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